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In this article we study the thermal response functions for two one-dimensional models, namely 
the Hubbard and spin-less fermion t-V models. By exactly diagonalizing finite sized systems we cal- 
culate dynamical electrical, thermoelectrical, and thermal conductivities via the Kubo formalism. 
The thermopower (Seebeck coefficient), Lorenz number, and dimensionless figure of merit are then 
constructed which are quantities of great interest to the physics community both theoretically and 
experimentally. We also geometrically frustrate these systems and destroy integrability by the inclu- 
sion of a second neighbor hop. These frustrated systems are shown to have enhanced thermopower 
and Lorenz number at intermediate and low temperatures. 
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I. INTRODUCTION 

Strongly correlated electron systems are currently at 
the forefront of physics. These systems provide some 
of the most interesting theoretical challenges and many 
experimental systems of fundamental and technologi- 
cal interest consist of strongly correlated constituent 
particles 1 ' 2 . It is at the convergence of the fundamen- 
tal and the technological that thermal response functions 
of strongly correlated systems take on particular impor- 
tance. This is certainly true in the high thermopower 
material sodium cobalt oxide 3 ' 4 ' 5 ' 6 as well as other tran- 
sition metal oxides such as the high T c superconduc- 
tors. Thermal behavior is also of interest in experimen- 
tal one-dimensional systems such as carbon nanotubes 7 , 
semiconductor nanowires 8 , and organic compounds 9 ' 10 , 
to name a few. 

One-dimensional systems are very interesting from the 
point of view of the exotic collective behavior of excita- 
tions due to reduced dimensionality. This is manifested 
most strikingly in the existence of a Luttinger liquid of 
bosonic excitations in generic ID systems at low temper- 
atures. Electrical and thermal transport in these systems 
is of great interest 11 . One-dimension also allows the ex- 
istence of a class of systems known as integrable systems, 
where there is an infinite family of commuting operators 
that commute with the Hamiltonian 12 ' 13 ' 14 . Recently, 
Shastry 15,16,17 has introduced a high frequency formalism 
for thermoelectrics (discussed below) that is particularly 
suited for strongly correlated electron systems. Hence, 
one-dimensional systems provide good systems in which 
to benchmark his new high frequency formalism. 

Geometrical frustration 18 has generated much inter- 
est in the physics community throughout past decades. 
Recently it has been shown that strong electron corre- 
lations in conjunction with the geometrical frustration 
induced by a two-dimensional triangular lattice are the 
keys to understanding the Curie- Weiss metallic phase in 
sodium cobalt oxide 6 . Geometrical and/or electrical frus- 
tration is also the key to the emergence of kinetic anti- 



ferromagnetism 19 , and to a description of quantum spin 
glasses 18 . 

In a one-dimensional system it is possible to induce 
frustration by considering kinetic energy hoppings fur- 
ther than nearest neighbors, i.e., second nearest neigh- 
bors. We demonstrate in this paper that frustration of 
this sort also has very interesting effects on transport 
properties in addition to the equilibrium properties al- 
luded to above. One such effect is enhancement of the 
thermopower compared to the unfrustrated system re- 
cently predicted by Shastry 6 ' 15 ' 16 ' 17 . Understanding this 
enhancement could lead to clues towards the creation, in 
material science laboratories, of custom made high ther- 
mopower materials with far reaching consequences to the 
scientific community. 

Theoretically, strongly correlated systems are notori- 
ously difficult due in part to the fact that perturba- 
tion theory is not applicable. For dynamical thermal 
response functions (using the Kubo 20 formalism) , knowl- 
edge of the full eigenspectrum is necessary, therefore, 
progress is made via exact diagonalization of finite sized 
systems. The models and systems we study are the Hub- 
bard and spin-less fermion t-V model on one-dimensional 
rings (periodic boundary conditions) of up to £ = 10 
sites for the Hubbard model and C — 16 sites for the 
spin-less fermion t-V model. Newer theoretical meth- 
ods are constantly being developed, such as the dynam- 
ical mean field theory 21,22 , finite temperature Lanczos 
(FTL) methods 23 , cluster perturbation theory 24 , etc.. In 
fact, the dynamical mean field theory approximation has 
considered some thermoelectric variables previously, al- 
though in a limited range of model parameter space, c.f. 
Ref. 25 and 26. It is important, as a first step, to ap- 
proach these systems with a rigorously exact method (ex- 
act diagonalization) to, if nothing else, provide a bench- 
mark for further approximate studies and methods. Of 
course, there are obvious shortcomings to exact diago- 
nalization, i.e., finite size effects due to the smallness of 
the system and the computational price is often quite ex- 
pensive. However, we find only the very low temperature 
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regimes (T < l|f|) troublesome in our studies which, inci- 
dentally, are the regimes that other approximations also 
have difficulty tackling. Interestingly, the high frequency 
formalism of Shastry, compared to dynamical approxima- 
tions, is significantly less expensive computationally and 
perhaps could lend itself even more profitably to new and 
existing approximate methods. 

We primarily use the Hubbard model for our studies 
and carry out the most extensive calculations on it. The 
t-V model is mainly used to supplement the study of 
the Hubbard model, partly due to the fact that some of 
the transport coefficients turn out to be identically zero 
in this model as discussed in Sec. IV B and partly be- 
cause the calculation of some of the operators are very 
involved, owing to the interaction being between neigh- 
boring sites and not on-site. The numerical calculation of 
the thermopower is the most tractable and is performed 
extensively. The calculations presented here are the first 
detailed calculations (considering the full range of model 
parameter space) in the literature for the thermopower, 
Lorenz number, and dimcnsionlcss figure of merit for 
strongly correlated models. 

The plan of this paper is as follows: in section II we in- 
troduce the models to be studied and provide an overview 
of our exact diagonalization procedure. Section III re- 
views the Kubo formalism for the thermoelectric conduc- 
tivities needed to calculate the physical observables of in- 
terest. Section IV, V, and VI present results for the ther- 
mopower, Lorenz number and thermoelectric figure of 
merit (FOM), respectively, for the models we study. We 
present our conclusions and a summary in section VII. 

II. MODELS 

The Hubbard model is described by a Hamiltonian 
with a kinetic energy term which allows electrons to hop 
between sites j and j + rj with probability t{rf) and an 
on-site electron repulsion potential energy governed by 
parameter U, i.e., 

c 2 

j=l r)=l a 
C 

+tr£(n iT -l/2)(n ii -l/2) ) (1) 
i=i 

where Cj a (cj a ) create (destroys) an electron with spin 

a — (|, |) at the lattice site j, hj a — $j a Cj a is the num- 
ber operator, and r\ — 1 and r\ = 2 for first and sec- 
ond nearest neighbor hoppings, respectively. The elec- 
tron operators obey the usual anti-commutation rules 
{cja,Cj, a ,} = 5 j: j'6 ar7 >. We assume t(l) = t, t{2) = t! , 
and t{if) — for all other rj. With t' = 0, this Hamilto- 
nian is particle-hole symmetric as can be seen from the 
local transformation dj a = (— l) J Cj cr . This transforma- 
tion leaves t and U invariant, while taking N a — > 1 — N a , 



where N a is the total number of electrons with spin 
a and the total number of electrons is N — 
Since it takes t' — ► — t', a non-zero t' breaks particle- 
hole symmetry. Even in the presence of t' ^ 0, this 
transformation is useful since it tells us that a quantity 
A{ V T ,N l ,t,t',U) = ±A(C -N h C-N h t,-t,U). Thus 
knowing the value of A up to half-filling for t' and —t' 
lets us construct the entire dependence on filling in both 
cases, provided we know whether A is odd or even un- 
der the particle-hole transformation. We only consider 
electron repulsion U > and assume periodic boundary 
conditions (one-dimensional rings). 

The spin-less fermion t-V models is governed by a 
Hamiltonian with a kinetic energy term similar to the 
Hubbard model (without the spin) but the potential en- 
ergy describes a nearest neighbor repulsion governed by 
a parameter V, i.e., 

c 2 

H = -^^(ryHcj^ + h.c.} 

j = \ r,= l 
C 

+vY,(n 3 -l/2)(n 3+1 -l/2) . (2) 

i=i 

Like the Hubbard model we assume here that i(l) = 
t, t(2) = t' and t(r?) = for all other 77 and V > 0. 
This model too has the same particle-hole symmetry but 
without the spin. Once again, t' breaks the particle-hole 
symmetry and a calculation up to half-filling (with t' > 
and t' < 0) lets us extrapolate to all values of filling. 

As with any exact diagonalization there exist compu- 
tational limitations due to the obvious difficulties in deal- 
ing with large matrices. In view of these limitations 27 the 
Hubbard model is computed for N =1,. . .,5 and 5,. . .,8 
electrons on C =10 and 8 site systems, respectively. 
These particular systems correspond to electron fillings 
(densities) of n = N/C =0.1, 0.2, 0.3, 0.4, 0.5, 0.625, 
0.75, 0.875, and 1. The t-V model allows a larger lattice 
size stemming from the lack of the spin degree of free- 
dom and we present data for L — 16 (and N = 1, . . . , 8). 
The large Hilbert spaces can be reduced by consider- 
ing systems with a constant z-component of spin (for 
the Hubbard model). The sector with the smallest z- 
component of spin has the largest Hilbert subspace di- 
mension and dominates the physics. Linear momentum 
is a good quantum number due to the translational in- 
variance of our systems and we implement this symmetry 
to further reduce the Hilbert space to more manageable 
proportions. However, we point out that the bottleneck 
of our calculations is not necessarily the Hamiltonian di- 
agonalization but the finite temperature averages of cer- 
tain operators and evaluation of the full Kubo formulas 
discussed below. 

The density of states of our models for t' = has two 
van Hove singularities at the band edges. If \t'\ > 0.25|i| 
another singularity is introduced which will invariably 
produce large changes in the thermoelectric properties. 
However, we are interested here in the more subtle ef- 
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fects arising from geometrical frustration and strong in- 
teractions, hence, by introducing a non-zero \t'\ < 0.25|f | 
the density of states maintains only the original two van 
Hove singularities. 

In order to display effects of electron correlations the 
interaction parameter U(V) must be at least a few times 
larger than the bandwidth which is equal to W = 4|t| 
for both models. We consider three values of interac- 
tion strength for each model corresponding roughly to 
the non-interacting case where the interaction parameter 
is equal to zero (U(V) = <C W), the weakly interact- 
ing case where the interaction is equal to roughly the 
bandwidth (U(V) ~ W), and the strongly interacting 
case where the interaction is a few times greater than 
the bandwidth (U(V) > W). 



III. DYNAMICAL THERMAL RESPONSE 
FUNCTIONS 

Conductivities arc computed via the Kubo linear re- 
sponse formalism 28 recently presented in detail by Pe- 
terson et al in Ref. 29 which closely followed the work 
of Shastry 15,16 ' 17 . In particular, we calculate the Kubo 
formulas for the electrical thermoelectrical 7(w), 

and the thermal k(uj) conductivities, respectively, which 
then allows us to calculate physical quantities of interest 
such as the thermopower S (or Seebeck coefficient), the 
Lorenz number L, and the dimcnsionlcss figure of merit 
ZT, given as 



S(oj,T) = 



(3) 
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FIG. 1: Probability density of states P(s) with energy dif- 
ference s = \si — Sj\. The solid black and dashed blue lines 
represent U = 4|£| and 12\t\, respectively, for density n = 0.3 
(a), n — 0.5 (b), and n — 0.9 (c). The position of the maxi- 
mum probability for the lower Hubbard band is insensitive to 
the value of U. 



L(cj,T) 



k(ui,T) 
Ta(u), T) 



and 



Z{uj,T)T = 



L(u,T) 



(4) 



(5) 



given here for completeness and ease of presentation. The 
second term in Eq. 4 is produced by the zero electric cur- 
rent constraint under a thermal gradient 30 . This term is 
usually small (especially at low temperatures) for met- 
als and semiconductors and often ignored. However, as 
will be shown, for strongly correlated lattice systems like 
those studied here, this term is an important actor espe- 
cially at high temperatures. 

However, one is usually interested in the DC (oj — > 0) 
limit of the dynamical conductivities. For finite sized 
systems this requires the introduction of a level broad- 
ening e which smoothens the divergences caused by the 
discrete nature of the eigenspectrum. Often e is taken 
to be equal to the mean energy level spacing of the sys- 
tem which is of order 0(\t\/C). However, for the Hub- 
bard model, as the interaction energy U is increased the 



mean energy level spacing begins to include upper Hub- 
bard bands. In Fig. l(a)-(c) the probability density of 
states P(s) with energy difference s — |e$ — £j\ (ek is the 
eigenenergy of state k) is plotted for three representative 
cases, namely, fillings n = 0.3, 0.5, and 0.9, for U = 4\t\ 
(weakly-coupled) and U = 12|i| (strongly-coupled). 

It is clear from Fig. l(a)-(c) that the most probable 
energy difference between states in the lower Hubbard 
band is relatively immune to changes in cither filling or 
interaction strength. However, the appearance of the up- 
per Hubbard bands in the strongly-coupled case yields 
an e that is very large (compared to the bandwidth) and 
strongly U dependent. The large U situation yields a 
large e which tends to mask real physical contributions 
to the current matrix elements coming from transitions 
between Hubbard bands. Therefore, in this work we 
choose e to be approximately equal to the mean energy 
level spacing in the lower Hubbard band for all cases, i.e., 
e = 0-8 1 £ | for the Hubbard model 31 . 

As discussed previously in detail in Ref. 29 another 
frequency limit, besides the DC limit, is the infinite fre- 
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quency limit which is defined as (for completeness) 
S*(T) = lim S(w,T) '''' ' 



T (f x 



(6) 



L*(T) = lim L(w,T) 



T2(f. 



[S*(T)] 2 , (7) 



and 



Z*(T)T= lim Z{u,T)T = 



L*(T) 



(8) 



with the operators f xx , $ xx , and Q xx defined in 15,16,17 . 
These quantities can be calculated as equilibrium expec- 
tation values of operators that, while non-trivial, are eas- 
ier to calculate and less time consuming numerically than 
the full dynamical quantities via the Kubo formulas. All 
of the interaction effects remain in these quantities but, 
importantly, the dynamics have been separated from the 
interactions. 



IV. THERMOPOWER 

The thermopower, or Seebeck coefficient, is defined as 
the ratio of the thermoelectrical to the electrical conduc- 
tivity and measures the electrical response of a material 
to a temperature gradient. Like the Hall coefficient it is 
often assumed to measure the sign of the charge carriers 
in a particular material. 

The thermopower S(u>, T) can be instructively 
rewritten 6,29 , by isolating the term containing the chem- 



ical potential fJ-(T), as 



S(cj,T) = S tr (u,T) + S M H(T) 



(9) 



The first term is due to electrical transport while the sec- 
ond is entropic in origin and is the familiar Mott-Heikes 32 
term S M h{T) = -0( T ) - M°))A?eT. These two terms 
both contribute to the thermopower in different and often 
competing ways. 

Strongly correlated systems are difficult to make the- 
oretical progress on and the first term in Eq. 9 is usu- 
ally the most intractable. For many systems this term is 
small (especially at low temperatures), compared to the 
second term, and fruitfully ignored. In this approxima- 
tion the thermopower is dominated by the Mott-Heikes 
(MH) term. The MH limit is described as the limit when 
ksT 3> This limit is achieved at either very high 
temperatures or at more modest temperatures in nar- 
row band systems 32 . The usefulness of the MH term is 
that at high temperatures the chemical potential is lin- 
ear in T leaving the MH limit constant. Even though 
the MH limit is constant it contains a non-trivial filing 
dependence that arises from the particular nature of the 
Hilbert space. Many of these MH limits were considered 
previously 32 and below we quote them for the spin-less 



fermion t-V model for finite and infinite V in Eq. 11. As 
will be shown below, the transport term (first term in 
Eq. 9) is very important for low to intermediate temper- 
atures and the MH term alone is clearly inadequate in 
this temperature regime. 

The thermopower is expected to vanish in the zero 
temperature limit. This is physically intuitive and has 
been shown theoretically for the one-dimensional Hub- 
bard model via a version of the Bethe ansatz solution 33 . 
In our formalism this T — vanishing is accomplished 
through a subtle balance of the transport and MH terms 
as T — > 0. For thermodynamically large systems this bal- 
ance is obtained. However, for finite systems, even in the 
non-interacting limit, the balance is not manifest result- 
ing in thermopower divergences as T — > due to finite 
size effects. 

This balance suggests a rewriting of the thermopower, 
given in Ref. 6,29, in such a way to ensure that the ther- 
mopower vanishes at T — by forcing both the transport 
(S tr (uJ,T)) and MH (S M h{T)) terms to independently 
vanish in that limit. This rewriting yields a frequency de- 
pendent transport term and the frequency independent 
MH term. 

Finite sized systems have a few more subtleties which 
we now describe. The chemical potential is commonly 
defined as fi{T) = dF N (T)/dN, where F N (T) is the 
Hclmholtz free energy for an N particle system. For our 
finite systems we will approximate this partial deriva- 
tive as n(T) = (F N+1 (T) - Fn_i(T))/2 for N > 2 or 
fj,{T) = F N+1 (T) - F N (T) for N = 1 since there is no 
system for N = 0. As discussed in Ref. 29 the ground 
state degeneracy of a system (if it exists) is discounted 
when calculating /i(T) so as to eliminate a leading order 
term linear in T that produces a non-zero Smh(T) as 
T — ► 0. Further, finite systems have discrete energy lev- 
els giving rise to an energy gap. The energy gap causes 
an exponential behavior in ^t(T) which is not a serious 
problem since it vanishes faster than T 2 . Both of these 
particulars, however, create a chemical potential which 
does not behave as T 2 at low T as expected for thermo- 
dynamically large systems. 

In the following figures the thermopower will be given 
in its "natural" units of kg/q e where q e = — |e| with 
|e| the value of the electron charge. For experimental 
comparison one simply replaces fcs/|e| = 86/iV/K. 



A. Hubbard model 

The MH term is the high temperature limit ksT ^> \t\ 
of S M h(T). As discussed elsewhere 29 ' 32 the MH limit of 
the finite U situation is essentially the uncorrelated band 
since the temperature is necessarily much larger than U. 
To understand the effects strong interactions play (large 
U) one must consider infinite U when calculating this 
limit. 

In the finite U case, the MH limit has a single zero 
crossing at half filling n = 1 where the thermopower is 
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seen to change sign and it diverges in the band-insulator 
limits (n = and n — 2). For the infinite U case 
there are two additional zero crossings and one addi- 
tional divergence. The MH limit still diverges in the 
band-insulator limits but now also diverges for the Mott 
insulator (n = 1). The two additional sign changes occur 
for fillings n = 2/3 and n — 4/3. The MH limit in the 
infinite U case for n > 1 is found from the MH limit for 
> n > 1 through particle-hole symmetry (n — > 2 — n 
and q e — > — q e ). There is no t' dependence for the MH 
limits, therefore, they lead to the conclusion (for t and 
if) that the Hubbard model should have interaction in- 
duced sign changes of the thermopower at n « 2/3 (and 
n rs 4/3 using particle-hole symmetry). 

First we consider results for the one-dimensional Hub- 
bard model with t' — 0. 



i. Hubbard with t = 

Here we consider the Hubbard model with t' = 
in the non-interacting (U = 0), the weakly correlated 
(U = 4|t|), and the strongly correlated (U = 12|t|) cases, 
respectively. 

Due to the particle-hole symmetry of the if — Hub- 
bard model /x(T) for N — C is 



1 , /^-i 
273 ln fe 



1 



In 



£-1 



17 
2~ 



(10) 



since an energy eigenstate for N electrons and 2C — N 
electrons are related through Sk(N) = £fc(2£ — N) — (C — 
N)U according to Lieb and Wu 34 . Hence, Smh(T) is 
identically zero for all T. The transport term St r (T) is 
also identically zero for all T due to particle-hole symme- 
try and hence we find that S(u>, T) = for n = 1 for all 
T and lo. This argument equally applies to the spin- less 
t-V model for t' = 0. 

In Fig. 2(a)-(c) we plot the infinite frequency limit of 
the thermopower S*(T) and the DC limit (uj — > 0) of 
the full Kubo thermopower S(u), T) versus temperature 
T and filling n for U = 0\t\ (a), 4|t| (b), and 12\t\ (c). 
Projected onto the the T = 5\t\ plane are the MH limits 
for finite and infinite U scenarios. 

For the non-interacting case in Fig 2(a) there is no fre- 
quency dependence as both the charge and heat current 
operators are diagonal. The thermopower grows mono- 
tonically and, essentially, linearly from zero at T = to 
the MH limit 35 near T — 5\t\. As the temperature is 
increased the slope of the linear regime is lessened and 
the thermopower approaches the MH limit more slowly. 
Some finite size effects are evident at the largest fillings 
n =6/8 and 7/8 calculated where some somewhat strange 
low temperature effects are observed. 

The effect of interactions shown in Fig. 2(b)-(c) is to 
generally reduce the thermopower for fillings greater than 
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FIG. 2: (color online) S*(T) (black line) as a function of filling 
n and temperature T. For (b) and (c) the orange dotted 
lines are the DC limit 5(0, T) of the full thermopower for 
comparison. Projected onto the T = 5[i| plane are the MH 
limits for both the finite (red line) and infinite U (blue dashed 
line) situations. 



n = 2/3. For fillings below n ~ 0.625 the DC and infi- 
nite frequency thermopower arc essentially identical. At 
fillings greater than n = 0.625 S(0, T) and S*(T) dis- 
play slight differences at low temperatures for the weakly 
correlated regime (U = 4|t|, Fig. 2(b)) and marked dif- 
ferences for the strongly interacting regime (U = 12|t|, 
Fig. 2(c)). 

The thermopower at half filling in all cases is pinned at 
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t7t=0, U=12W 




(c) 



FIG. 3: (color online) Smh(T) (black line) as a function of 
filling n and temperature T. For (b) and (c) the orange dotted 
lines are the infinite frequency limit S*(T) for comparison. 
Projected onto the T — 5\t\ plane are the MH limits for both 
the finite (red line) and infinite U (blue dotted line) cases. 



zero as discussed. For n = 0.875 the effect of interactions 
is to reduce the thermopower (Fig 2(b)), changing its 
sign (Fig. 2(c)), for low to intermediate temperatures. 
Eventually the entropy dominates and the thermopower 
begins to climb towards its entropy determined MH limit. 
It should be noted that even for the highest temperature 
T = 5\t\ the thermopower remains negative for n = 0.875 
and U = I2\t\. 

In Fig. 3(a)-(c) we plot the MH term of the ther- 



mopower (Smh(T)) versus temperature and filling for 
U = 0, 4|i|, and 12|£|, respectively. Also shown is 
S*(T) for comparison. Interestingly, the transport term 
increases the low to mid temperature thermopower for 
all fillings in the non-interacting and weakly interacting 
cases. For the strongly interacting case in Fig. 3(c) the 
transport term begins to decrease the thermopower at 
low temperatures and high fillings (n > 0.75). 

The results in this section can be directly compared 
with Ref. 36 which used the FTL method which allowed 
rings of up to C — 14 sites. Their results were for a much 
lower window of temperatures 0.2 \t\ < T < 2\t\ and show 
a greater thermopower suppression for more modest val- 
ues of U than in the present work. The thermopower 
in the low temperature regime is very sensitive to any 
changes in either the S tr (T) or Smh{T). Considering 
that Ref. 36 used a larger system but an approximation 
it is unclear as to the nature of the discrepancy between 
those results and the present study. However, the quali- 
tative behavior of both calculations is consistent. 

Next the full frequency and temperature dependence of 
the thermopower is investigated. Fig. 4(a)-(h) plots the 
difference in the full frequency dependent thermopower 
minus the infinite frequency limit, i.e., S(lu,T) — S*(T). 
For the weak coupling (Fig. 4 (a) -(d)) case there is lit- 
tle frequency dependence except for when w ~ U at low 
to intermediate temperatures. In what is surely a con- 
sequence of the finite sized system there appears to be 
an even/odd effect in that for even numbers of electrons 
(Fig 4(c) and not shown) S(u),T) - S*(T) is positive for 
u! < U and negative for u> > U while for odd numbers 
of electrons (Fig 4(a)-(b)-(d)) the opposite effect is ob- 
served. Of course, for half filling there is no frequency 
dependence (not shown). 

The same qualitative behavior is found for the strongly 
correlated regime in Fig. 4(e)-(h). Again the even/odd 
effect is obtained (cf. Fig. 4(g) vs. Fig. 4(e), (f), and (h)) 
and the frequency dependence is generally weak except 
for u> ~ U . The main difference between the U = I2\t\ 
and the U = A\t\ cases is that the frequency dependence 
at lu ~ U in the former is larger. It should be kept 
in mind that the largest frequency dependences occur for 
the lowest temperatures and the thermopower for a finite 
sized system in the low temperature regime is suscepti- 
ble to finite sized effects that could possibly not survive 
in the thermodynamic limit. Therefore, it is safe to as- 
sume that our results are qualitatively representative of 
a thermodynamically large Hubbard model but perhaps 
not quantitatively accurate. 

Similar qualitative frequency dependence was seen in 
Ref. 36 for a low temperature slice (T = 0.5|t|) and high 
fillings (n > 0.7). The difference in that work and the 
present work is that the frequency dependence occurred 
nearer to w ~ U/2 in the former where in the present 
result the largest frequency dependence occurs near lu > 
U. 

Unless one is concerned with extremely low temper- 
atures and frequencies similar in magnitude to the in- 
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(S(0))-S 
0.4 
0.3 
0.2 
0.1 

o 

-0.1 
-0.2 





(c) 



(d) 



n=0.3, Wt=0, U=12lll 



n=0.5, t7t=0, U=12W 



(S(o»-S )[k B /qJ 

0.3 
0.2 
0.1 



(S((0)-S )[k B /q e ] 





(e) 



(f) 



n=0.75, l'/t=0, U=12lll 



n=0.875. l'/l=0, U=12ltl 



(S(0))-S )[k B /qJ 





(g) 



(h) 



FIG. 4: (color online) S(cu,T) — S*(T) as a function of frequency u) and temperature T for the weak coupling case U = 4\t\, 
t' — 0, and fillings (a) n = 0.3, (b) n = 0.5, (c) n = 0.75, and (d) n — 0.875. Panels (e)-(h) represent the strong coupling 
case U = 12|t| for the same fillings as the weak coupling case. The main frequency dependence occurs near uj ~ U for low to 
intermediate temperatures. The largest frequency dependence near T = is likely a finite size effect. There is no frequency 
dependence for the half filled case n — 1 (not shown). 
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(c) (d) 

FIG. 5: (color online) S*(T) (black line) as a function of filling n and temperature T for (a) (7 = 0, (b) U = 4\t\, and (c) 
U = 12|t| and t'/t = -02.. The full S(u,T) (not shown 37 ) compares similarly to S*(T) for t' ^ as for t' = 0, c.f. Fig. 4. The 
red curve is for t' — to facilitate an easy comparison. As an example of the other sign of t' , namely, t'/t = 0.2, (d) shows the 
expected reduction in S*(T) for U = 4\t\ (similar results obtained for (7 = and U = 12\t\ are not shown). The MH limits are 
projected unto the T = 5\t\ plane for finite U (red line) and infinite U (blue dotted line). 



teraction strength U the infinite frequency thermopower 
S*(T) is a good representative of the full thermopower 
S(w,T). Recently, a similar calculation on the strongly 
correlated t-J model for a two-dimensional triangular lat- 
tice was carried out in Ref. 29 where it was found that 
the frequency dependence was much weaker justifying the 
use of S*{T) in place of S{w, T). 



2. Hubbard with t' =£ 

Here we investigate the effect of frustration on the ther- 
mopower by considering the case of a non-zero second 
neighbor hopping amplitude \t'\/t — 0.2. 

Recently, Shastry 15,16 ' 17 predicted a low to intermedi- 
ate temperature enhancement of the thermopower via a 
high temperature expansion of the high frequency limit 
S*(T) for the geometrically frustrated two-dimensional 
triangular lattice t-J model. In that case, changing 
the sign of the hopping was found to enhance the ther- 
mopower at intermediate temperatures. Haerter et. al. 6 
and Ref. 29 more thoroughly investigated that particular 
case. 



In one-dimension it is similarly expected that an en- 
hancement of the thermopower will occur when a second- 
neighbor hop is added to the kinetic energy which frus- 
trates the lattice and destroys integrability. 

Fig. 5(a)-(d) show S*(T) versus temperature and filling 
for values of t'/t = -0.2 (Fig. 5(a)-(c)) and t'/t = 0.2 
(Fig. 5(d)). For this case we do not plot the full frequency 
dependence but remark that it is similar qualitatively and 
quantitatively to the t' = case . 

When t > and t' < the thermopower is enhanced at 
low temperatures. The enhancement is seen to arise from 
almost purely the transport term of the thermopower 
(first term in Eq. 9). Fig. 6 shows S*(T) and the MH 
term S MH (T) for t'/t = -0.2 and U = 4|i| as a rep- 
resentative example. Similar to the results displayed in 
Fig. 3 the low temperature thermopower is dominated by 
the transport term and in the case of t'/t = —0.2 that 
domination is even more pronounced as it produces an 
enhancement peak. 

For small fillings the enhancement has weak U de- 
pendence although the non-interacting case seems to be 
tending to diverge at very low temperatures (n = 0.1 
and n — 0.2 in Fig. 5(a)) that is most certainly a finite 
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f/t=-0.2, U=4ltl 




FIG. 6: (color online) Smh(T) (black line) as a function of 
filling n and temperature T. The orange dotted lines are the 
infinite frequency limit S* (T) for comparison. Projected onto 
the T — 5\t\ plane are the MH limits for both the finite (red 
line) and infinite U (blue dotted line) cases. 



size artifact. At fillings above n = 0.5 there are inter- 
action effects which are visible at low temperatures. For 
the frustrated case discussed here the thermopower still 
mostly achieves its MH limit by T — 5|i| as expected. 

For the situation when t and t' are both positive there 
is a suppression of the thermopower at low temperatures, 
c.f. Fig. 5(d), hence the opposite effect. Only the weakly 
correlated U = 4|t| case is plotted to illustrate this. 

For both signs of t' at half filling (n = 1) the ther- 
mopower is no longer identically zero since the addition 
of a non-zero t' destroys the particle-hole symmetry, how- 
ever, the thermopower remains quite small at this den- 
sity. 



B. t-V model 

Most of the considerations of the thermopower for the 
Hubbard model also apply to the t-V model. In partic- 
ular this model too has particle hole symmetry at half- 
filling (n = 0.5 instead of n = 1 as for the Hubbard 
model, since the spin degree of freedom is absent) and 
thus the thermopower is identically zero at all tempera- 
tures. The thermopower is once again divergent in the 
band-insulator limits (n = and n = 1) and also in the 
vicinity of the Mott insulator (n = 0.5) in the presence 
of strong correlations (large V). In the high tempera- 
ture limit (ksT 3> t), the thermopower is pinned by the 
entropy as in the Hubbard model. The MH limit can 
be calculated in a straightforward manner in the V = 
and V — oo regimes from the expressions already de- 
rived for the extended Hubbard model (on-site U and 
nearest-neighbor V) 32 . The t-V model corresponds to 
setting U = oo and removing the spin degeneracy. We 



f/t=0, V=0ltl 



(a) 



(b) 




(c) 

FIG. 7: (color online) S*(T) (black line) as a function of filling 
n and temperature T. For (b) and (c) the orange dotted 
lines are the DC limit 5(0, T) of the full thermopower for 
comparison. Projected onto the T = 5[i| plane are the MH 
limits for both the finite (red line) and infinite V (blue dotted 
lines) situations. 



thus obtain 



lim S M h(T) = 




where Eq. ll(i) corresponds to finite V and < n < 1 
while Eq. 11 (ii) corresponds to infinite V and < n < 
0.5. As mentioned previously, the expressions in Eq. 11 
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(c) (d) 

FIG. 8: (color online) S*(T) as a function of filling n and temperature T for (a) V = 0, (b) V = 8|t|, and (c) V = 12\t\. The 
full S(u>, T) compares similarly to S*(T) as for the t' = case, c.f. Fig. 4. The red curve is S*(T) for t' = to facilitate an easy 
comparison. As an example of the other sign of t' , namely, t' /t = 0.2, (d) shows the expected reduction in S*(T) for V = 8\t\. 
The MH limits are projected unto the T = 5\t\ plane for finite V (red line) and infinite V (blue dotted line). 



were first considered in Ref. 32 with the the second line 
(the infinite V limit) being the third equation of Table I 
in Ref. 32 with the spin degeneracy removed. 

An important distinction between the Hubbard model 
and the t-V model is that the energy current operator 
commutes with the Hamiltonian in the latter. From the 
Kubo formula, this implies that in Eqn. 9, S tr (uj,T) = 
for all uj j= and S(w, T) = S MH (T) for all lo ^ 0. Since 
S*(T) = S(uj -> oo, T), this leads to S*(T) = S M h(T) 
at all temperatures and filling in this model. From 
Eqn. 6, we have that = —fJ-(t xx ), a fact that has 

been verified numerically. A further consequence is that 
k(lu,T)/T<7(lu,T) = S(lu,T) 2 for all u ^ and conse- 
quently L(u > 0,T) = implying Z(u > 0,T)T = oo 
in this system. Physically, the energy current commuting 
with the Hamiltonian means that the system is unable to 
transport any heat at finite frequency in the presence of 
a temperature or potential gradient without transport- 
ing charge. In zero current conditions where there is no 
charge transport, there is no heat transport as well and 
the thermal conductivity is zero. This causes the Lorenz 
number to be zero and the FOM to be infinite. All the 
above considerations rely on the fact that the energy cur- 
rent commutes with the Hamiltonian. This is no longer 



true with t' ^ and all the quantities mentioned above 
will have have non-zero values at u> 0. 



1. t-V with t' = 

The behavior of the thermopower for the t-V model is 
quite similar to that of the Hubbard model even though 
the thermopower has no transport term for t' = and 
is, therefore, simply Smh(T). Plotted in Fig. 7(a)-(c) is 
the thermopower S*(T) for three values of the interac- 
tion strength, V = 0, V = 8\t\, and V = 12\t\. At low 
densities the thermopower very rapidly rises to nearly its 
full MH limit by T w 2\t\. The initial (mostly) linear 
slope is reduced as the filling is increased. The effects of 
interactions are not readily seen until the somewhat large 
filling of n — 0.4375 where the thermopower is markedly 
reduced for all temperatures shown. Eventually, as the 
temperature is raised, the interaction effects are washed 
out as the thermopower approaches its MH limit. 

Compared to the Hubbard model, however, the inter- 
action effects appear weaker (perhaps due to the absence 
of the transport term) and effect only the highest fillings 
studied (other than the half filled case n = 0.5 which is 
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pinned at zero due to particle-hole symmetry). 

2. t-V with t' ^ 

Adding a second-neighbor t' hopping term has an ef- 
fect very similar to the one in the Hubbard model again 
destroying the integrability and the particle-hole sym- 
metry of the model. A result of the latter is that the 
thermopower is no longer identically zero at half-filling. 
As in the Hubbard model, a more interesting aspect of in- 
troducing a second-neighbor hopping is that it produces 
frustration depending on its sign. 

The thermopower is plotted as a function of tempera- 
ture and filling for V = 0, V = 8|t|, and V = I2\t\ for 
the case of t'/f = -0.2 in Fig 8(a)-(c), and for V = 8\t\ 
for the case of t'/t = 0.2 in Fig 8(d). 

A positive sign of t' reduces the value of S* (T) com- 
pared to t' — 0, while a negative sign enhances it. S*(T) 
starts out being zero at T = and approaches the Mott- 
Heikes limit as T — > oo independent of the value of t! . 
Since S* (T) for t! = is a monotonically increasing func- 
tion of T, it stands to reason that S* (T) reaches a maxi- 
mum at some value of T for t' < and decreases towards 
the Mott-Heikcs limit. This is indeed what is seen in our 
calculations, as demonstrated in Fig. 8, similar to the 
situation of the Hubbard model. This is very interesting 
since it affords the possibility of thermopower enhance- 
ment through frustration. 

Even though we do not explicitly calculate S(u>, T) for 
the t-V model when t' ^ 0, we comment that the in- 
troduction of t' causes the energy current operator not 
to commute with Hamiltonian. Consequently S(uj, T) is 
now no longer independent of u> and approaches S* (T) 
as lo — > oo like in the Hubbard model. Presumably, the 
full frequency dependence is similar to that seen in the 
Hubbard model but has not been calculated explicitly. 

Although not shown for the t-V model, we point out 
that similar to the Hubbard model, the thermopower 
enhancement peak arises from the transport term of 
the thermopower almost exclusively as discussed in 
Sec. IV A 2 and shown in Fig. 6. 

V. LORENZ NUMBER 

The Lorenz number (Eq. 4) is an important quantity as 
it measures the ratio of the thermal to electrical conduc- 
tivity. Further, it is a key ingredient to the FOM (con- 
taining both the thermopower and Lorenz number) which 
measures the efficiency of a thermoelectric material. Note 
that there is usually a contribution to the Lorenz number 
coming from lattice vibrations (phonons). In this work, 
however, we only consider the electronic contributions to 
the Lorenz number. 

As noted in Ref. 29 the chemical potential is absent 
from Eq. 4 and the Lorenz number can be understood 
completely within the canonical ensemble. Evidently, it 



t7t=0, U=0ltl 



(a) 



(b) 




(c) 

FIG. 9: (color online) L*(T)/Lo as a function of filing n and 
temperature T (black line). The orange dotted line is the 
DC limit L(0,T) for comparison. Note that for U = (a) 
there is no frequency dependence of L(u),T) and L*(T) = 
L(Q,T). The full L(u),T) (not shown) compares similarly to 
the frequency dependence of the thermopower for U — 4|t| 
and 12|t| in Fig. 4. 



is determined through electron transport alone. At zero 
temperature it is well known that the Lorenz number is 
equal to the constant L = (nkB/\/3q e ) 2 which is just 
the Wiedemann- Franz law 30 ' 38 . 

As discussed in Ref. 29, for non-interacting systems, 
the limit linty—o L(lu, T) = L comes from two effects. 
Similar to the thermopower, we attempt to limit the finite 
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(c) (d) 

FIG. 10: (color online) L*(T)/L as a function of filling n and temperature T for (a) U = 0, (b) U = 4|t|, and (c) U = 12\t\. 
The full L(u!,T) compares similarly to frequency dependence of the thermopower. The orange dotted curve is the t' = case 
to facilitate an easy comparison, (d) is an example of the other sign of t', namely, t' /t = 0.2. 



sized system induced divergences of the Lorenz number 
by forcing the subtle balances necessary to ensure a finite 
Lorenz number as T — > 0. The exact method used here 
is the same as in Ref. 29 and not repeated. Ultimately, 
we are able to control the divergences to a large degree, 
however, we are unable to answer the very interesting 
question of whether the value of Lq is a universal constant 
independent of electron interactions. 

Below we present L* (T) as a function of temperature 
and filling. The frequency dependence of L(u),T) (not 
shown here) is comparable to S(oj, T) in that it is gener- 
ally weak with a feature near u> = U . 



the Lorenz number starts at L and quickly decreases 
as a function of temperature similar to what is shown 
here. In the weakly coupled regime (Fig. 9(b)) L*(T) 
is very similar to the non-interacting case for fillings be- 
low approximately n = 0.4. For n > 0.5, however, it 
is enhanced at low temperatures. The DC limit is quite 
similar to the infinite frequency limit showing the useful- 
ness and accuracy of L*(T) in the weakly coupled regime. 
When the interactions are strong (Fig. 9(c)) there is a 
much stronger enhancement of L*(T) for n > 0.3 and 
the enhancement persists to higher temperatures. The 
DC limit in this case is not as similar to L*(T) as it is 
for the weakly coupled regime. 



A. Hubbard with t' = 

Fig. 9(a)-(c) shows the "normalized" high frequency 
Lorenz number L*(T)/Lq as a function of temperature 
and filling for the non-interacting, the weakly coupled, 
and the strongly coupled situations for the case of 1/ = Q. 
In Fig. 9(b)-(c) the DC limit of the full frequency depen- 
dent Lorenz number is also plotted, i.e., L(0,T) / Lq. 

For the non-interacting case (Fig. 9(a)) L*(T) is sup- 
pressed at low temperatures as the filling increases to- 
wards half filling. For a thermodynamically large system 



B. Hubbard with i'/0 

For non-zero t' we plot L* (T) /Lq as a function of filling 
and temperature along with L*(T)/L for the t' = to 
ease comparison in Fig. 10(a)-(d). 

For the same sign of t' < that produced the ther- 
mopower enhancement the Lorenz number in Fig. 10(a)- 
(c) is also enhanced compared to the t' — situation at 
low temperatures. This enhancement is also more pro- 
nounced for low than for high fillings and increases with 
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increasing interaction strength U, especially close to half 
filling. 

When t' > in the weakly coupled regime (Fig. 10(d)) 
the Lorenz number is very slightly suppressed compared 
to the t' = situation as one would expect from the 
t' < results and the previous investigation of the ther- 
mopower. 

Generally, for both t' = Q and t' ^ the Lorenz num- 
ber is appreciably below Lq for all temperatures above 
approximately T ~ 1.5|t|. 




VI. FIGURE OF MERIT 



(a) 



The figure of merit is the number that is most impor- 
tant when it comes to technological applications regard- 
ing thermoelectrics, as mentioned above, being a measure 
of thermoclcctrical efficiency. In our calculations we are 
ignoring the lattice contribution to the Lorenz number so 
our FOM calculated is that due to only electronic contri- 
butions. 

Using the "tricks" to handle finite size effects for S* (T) 
and L*(T), c.f. Ref. 29, we calculate the high frequency 
expansion of the FOM (Z*(T)T) given in Eq. 8. 

A large FOM can arise in essentially two ways. One 
way is for the thermopower, which is squared in the nu- 
merator, to be large. The other is for the Lorenz number 
in the denominator to be small. The Lorenz number 
has been shown in Sec V to tend to zero as the temper- 
ature tends to infinity. The thermopower, on the other 
hand, reaches a constant, and finite, MH limit as T — > oo. 
Hence, the electronic contribution to the FOM will even- 
tually grow to infinity as the temperature is increased 
without bound. 

In the following subsections we plot Z(uj,T)T as a 
function of T and filling n. 



(b) 




t'/t=0, U=12ltl 




Hubbard with t' = 



(c) 



Fig. ll(a)-(c) shows the FOM for the t' = Hubbard 
model for the non-interacting case, the weakly coupled 
case, and the strongly coupled case. Also plotted is the 
DC limit of the full frequency dependent FOM. In the 
non-interacting case (Fig. 11(a)) the FOM grows appar- 
ently quadratically in temperature with a coefficient that 
decreases inversely with filling. This behavior is seen in 
the weakly coupled and strongly coupled cases as well, 
c.f. Fig. ll(b)-(c), although the FOM is decreased more 
for lower fillings. 

The full frequency behavior of the FOM (not shown) 
is very similar to the high frequency limit, evidently 
because the differences in the two limits for the ther- 
mopower and Lorenz number largely cancel one another 
out. 



FIG. 11: (color online) Z*(T)T as a function of filling n and 
temperature T as the solid black line for U = (a), U = i\t\ 
(b), and U = 12|t| (c). The orange dotted curve is the DC 
limit of the full Z(lo, T)T for comparison, i.e., Z(0, T)T. Note 
that for U — (a) there is no frequency dependence. 



B. Hubbard with t' ^ 

In Fig. 12(a)-(c) and Fig. 12(d) we investigate the ef- 
fects of frustration on the FOM for the t' < and t' > 
cases, respectively. 

For the situation that produced a low temperature 
enhancement in the thermopower and Lorenz number 
(i' < 0) the FOM is suppressed compared to t' — except 
for the non-interacting case where they are very similar. 
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t7t=-0.2, U=Oltl 



f/t=-0.2, U=4ltl 





(a) 



(b) 



t'/t=-0.2, U=12ltl 



f/t=0.2, U=4ltl 





(c) 



(d) 



FIG. 12: (color online) Z*(T)T as a function of filling n and temperature T for (a) U = 0, (b) U = 4|t|, and (c) [7 = 12|t|. The 
full frequency dependent Z(u>, T)T compares similarly to the frequency dependence of the thermopower. The orange dotted 
curve is for t' = to facilitate an easy comparison. For an example of the other sign of t' we plot the situation t' jt = 0.2 in 
(d). 



The larger interaction strength has the effect of further 
suppressing the FOM especially for low fillings. 

Alternatively, the case of t' > that produced suppres- 
sion in both S* (T) and L* (T) has a very slight enhance- 
ment in the FOM shown in Fig. 12(d). This nicely il- 
lustrates the complicated way in which the thermopower 
and Lorenz number combine to produce the FOM. 



VII. CONCLUSION 

We have computed thermoelectrical properties of the 
Hubbard model and the spin-less fermion t-V model on 
one-dimensional rings investigating, in particular, the 
thermopower (Hubbard and t-V), Lorenz number, and 
figure of merit (Hubbard only). Our calculations are 
the first detailed calculations of these thermoelectric vari- 
ables in the literature for strongly correlated models. By 
adding a second neighbor hopping term with amplitude 
t' , both positive and negative, we were able to destroy the 
integrability of these models and induce frustration. For 
t! < the Hubbard and t-V model displayed an enhanced 
thermopower at low to intermediate temperatures. For 



the Hubbard model the Lorenz number was also found 
to have low temperature enhancements for t' < 0. How- 
ever, the FOM did not produce the same enhancement 
for t' < but instead a suppression for non-zero inter- 
action strength U. Although, the FOM was modestly 
enhanced by the opposite hopping t' > at low fillings. 

For the Hubbard model, the thermopower had a gener- 
ally weak frequency dependence other than a sometimes 
large feature near u> ~ U. This behavior was also ob- 
tained, but not shown here, for the Lorenz number and 
FOM. This has the consequence that the high frequency 
versions of the thermopower, Lorenz number, and FOM 
recently proposed by Shastry 15 ' 16,17 provide a good ap- 
proximation to the full dynamical quantities for most val- 
ues of the system parameters. 
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